Structure of vortex liquid phase in irradiated Bi 2 Sr 2 CaCu 2 08 -5 crystals 



N. Morozov a , M.P. Maley a , L.N. Bulaevskii a , V. Thorsm0lle a 
A.E. Koshelev 6 , A. Petrean fc and W.K. Kwok & 
a Superconductivity Technology Center, Los Alamos National Laboratory, Los Alamos, NM 87545 
b Materials Science Division, Argonne National Laboratory, Argonne, Illinois 60439 

(February 1, 2008) 

The c-axis resistivity in irradiated and in pristine Bi2Sr2CaCu20g_5 crystals is measured as a 
function of the in-plane magnetic field component at fixed out-of-plane component B± in the vortex 
liquid phase at T — 67 K. From this data we extract the dependence of the phase difference correla- 
tion length inside layers on B± and estimate the average length of the vortex line segments confined 
inside columnar defects as a function of the filling factor / = B±/B$. The maximum length, about 
15 interlayer distances, is reached near / w 0.35. 
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The structure of the vortex liquid in highly anisotropic 
layered superconductors with columnar defects (CDs) 
produced by heavy ion irradiation is one of the most 
intriguing questions in the current study of the vortex 
state in high temperature superconductors. For the most 
anisotropic P^S^CaC^Os-a (Bi-2212) superconductor 
without strong disorder, neutron scattering and Joseph- 
son plasma resonance (JPR) data provide evidence in 
favor of a pancake liquid with very weak correlations of 
pancakes in different layers When CDs are intro- 

duced into these crystals a large decrease in the reversible 
magnetization is observed B, indicating that pancakes 
are predominantly situated on the CDs, even in the liquid 
state. Recent studies of JPR and c-axis transport 
j^U in irradiated Bi-2212 crystals reveal enhancement of 
c-axis correlations in some interval of out-of-plane mag- 
netic fields, B±, below the matching field, at tem- 
peratures T ~ 70 K. In other words, pancakes at these 
fields and temperatures appear to form aligned segments 
of vortex lines inside CDs, while outside of this region 
they are better described as a liquid of weakly c-axis- 
correlated pancakes, as in pristine crystals. Both c-axis 
transport and JPR involve the flow of Josephson cur- 
rents and are thus sensitive to the misalignment of pan- 
cake vortices between adjacent planes. CDs aligned with 
the c-axis will promote intcrplanc pancake alignment in a 
region of temperature where pancakes remain largely lo- 
calized on CDs. Even in the liquid state, where pancakes 
are mobile, filling of the available CD sites should lead 
to enhanced c-axis correlation from statistical consider- 
ations alone. In previous c-axis transport measurements 
0, a dip in the magnetoresistance, p(B±), of a Bi-2212 
crystal was observed at a magnetic field corresponding to 
a filling factor of CDs, / = B ± /B$ w 1/3. 

So far there has been no technique available for pro- 
viding quantitative information on the degree of c-axis 
correlation or on the length of pancake line segments. In 
addition, it is unclear whether vortex interactions play 
a significant role in enhancing correlation in the liquid 
state. Previously, a method was proposed, but not yet 



realized, to extract the pancake density correlation func- 
tion using data for JPR frequency as a function of Bu 
at fixed B± In this Letter we determine the aver- 
age length of vortex line segments as a function of B±, 
i.e. of the filling factor /, from transport measurements 
on single crystals of Bi-2212 with and without CDs. For 
this we have developed a method for extracting the phase 
difference correlation function along the layer from mea- 
surements of the c-axis conductivity, <r c , as a function 
of the magnetic field component parallel to the layers, 
Bu, at fixed B±. The component B± establishes the vor- 
tex state to be studied, while the component Bu serves 
as a probe of this state, as described below. Knowing 
the phase difference correlation length, we estimate the 
pancake density correlation length along the c axis. 

In Josephson coupled superconductors in the presence 
of a c-axis current, voltage is induced by slips of the phase 
difference between layers, (p n ,n+i ( r i *) ; as described by 
the Josephson relation V^ in +i = (h/2e)ip n _ n+ i. Here n 
labels layers, r = x, y are coordinates in the ab plane, 
and t denotes the time. The c-axis conductivity in the 
vortex liquid state is determined by the Kubo formula 

a c (B ± ,B ll ) = (sJ$/2T) J alt J dvG{v,t), (1) 

G(r,t) = 2{sm<£„ in+ x(Q)0)sin<£n,n+i( r >*)) 

Rj (C0s[(p ni n + l(r,t) - (0,0)]), (2) 

where Jq is the Josephson critical current and s is the 
interlayer distance. Time variations of the phase differ- 
ence are caused by mobile pancakes fl0|| induced by B± 
and by mobile Josephson vortices induced by the parallel 
component Bu . In the lowest order in Josephson coupling 
we split [ip n ,n+i (0,0) — tpn,n+i ( r , t)] m t° t ne contribution 
induced by pancakes and that caused by the unscreened 
parallel component Bu . Assuming that Bu is along the x 
axis, we obtain a simple expression for the contribution 
of the parallel component to the phase difference: 

<£„,„ + i(0,0) -<p„ jn+ i(r,i) re 

[V?n,n+i(0,0) - tp n<n+1 (r,t)] Bll =o ~ 27rsB||j//$ - (3) 
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Inserting this expression into Eqs. (|l|) and ([2|) we obtain 
<t c (Bx,B\\) = (ttsJ^/T) J drrG(r,B ± )J (aB lir ), (4) 
where Jq(x) is the Bessel function, a = 2tts/^q and 

/>oo 

G(r,B ± ) = / dtG(r,t,B ± ). (5) 
Jo 

The function G(r, t) describes the dynamics of the phase 
difference caused by mobile pancakes. If g(B±,B\\) = 
&c(B±, B\\)/a c (Bj_,0) is known in the vortex liquid, the 
correlation function G(r) = G(r)jg$l/4TrsTa c (B_i_,0) 
may be found using the inverse Fourier-Bessel transform: 

G(r,B ± ) = j dB^g{B h B ± )J Q {aB\\r). (6) 

The correlation lengths R and R\ of this function are 
denned by the relations: 

Note that R\ is related to the coefficient of the B 2 term 
in the expansion of g(B±, By) in B?,: 

g(B ± ,B {l ) « 1 - [Trsi?!^)/^] 2 ^. (8) 

Thus R\(B±) can be obtained independently from data 
for cr c (B±, Bn) at small Bii. 

Such a procedure to obtain the function G(r,B±_) is 
valid if the vortex state depends weakly on the Josephson 
coupling and hence on the probe field Sii (which affects 
Josephson coupling). Then this method is nondestruc- 
tive. Let us check first under what conditions we can 
neglect the effect of Josephson coupling on the equilib- 
rium vortex state. The energy of Josephson coupling in 
the correlated area irR 2 is ttEqR 2 / 'A 2 , which should be 
much smaller than the temperature T to be treated as a 
perturbation §. Here£ (T) = / \ 2 ab {T) . For Bi- 
2212, with 7 « 300 and A afc (0) « 2000 A, the Josephson 
coupling in the correlated area is w 0.2T at the maxi- 
mum value of R found below and at T > 60 K. Thus the 
effect of Josephson coupling, and hence Sii , on the equi- 
librium vortex state may be neglected. For dynamical 
parameters, such as cr c , higher order terms in Josephson 
coupling describing dynamic screening of B^ omitted in 
Eq. (3), may be important. This will be discussed below. 

We anticipate that in the pancake vortex liquid state 
in pristine crystals the characteristic lengths R and R\ 
of the correlation function G(r) are of order of the inter- 
vortex distance, a = ($o/B±) 1 / 2 , because each pancake 
here is mobile and induces phase slippage as described in 
Ref. pj. In irradiated crystals with CDs we anticipate 
much bigger R and i?i if pancakes form long segments 



of lines inside CDs. Then only ends of segments con- 
tribute significantly to the phase difference tf n ,n+i ( r ) t) 
and lead to suppression of Josephson coupling and <7 C , 
while the effect of pancakes in neighboring layers inside 
the segments is much smaller. 

For our experiments, high quality Bi-2212 crystals 
(T c ~ 85 K) of about 1 x 1.5 x 0.02 mm 3 were used. 
The irradiation by 1.2 GeV U 238 -ions was performed on 
the ATLAS accelerator (ANL). According to TRIM cal- 
culations these ions produce in Bi-2212 continuous amor- 
phous tracks with diameter 4-8 nm and length 25-30 /im. 
Below, we present the results for the samples irradiated 
with a density of CDs corresponding to the matching field 
Bq, = 2 T and for a reference pristine sample. 

Our measurements of c-axis conductivity were carried 
out in a cryostat with two superconducting magnets pro- 
viding magnetic fields in orthogonal directions. The mag- 
nets are controlled independently and provide fields up to 
8 T in one direction and up to 1.5 T in the other direction. 
Samples can be oriented along the axis of either magnet. 
Misalignment of the crystal c axis with respect to the per- 
pendicular component B± was detected by the asymme- 
try of p c {B\\ ) with respect to the sign of Bu . We adjusted 
the direction of the field components, providing asymme- 
try below 5% at Bu > 2 T and below 10% at lower fields. 
The normalized conductivity <?(-Bii) was calculated using 
the average resistivity p c (Bn) = [p c (B\\) + p c (-B\\)]/2. 
Two silver contact pads were deposited on both sides of 
the sample using a mechanical shadow mask. The mask 
provided a clean surface rim ~ 0.1 mm from the sample 
edges and 25 /jm separation between current and poten- 
tial pads. The area of the contact was « 0.75 mm 2 for 
current and w 0.05 mm 2 for potential terminals. The 
resistance of the current pads at room temperature was 
w 2 fl. A current of 1 mA, driven through the sam- 
ple, provides an ohmic I-V regime. For very anisotropic 
Bi-2212, in our range of the magnetic field and tempera- 
tures, in-plane conductivity, a a b ~ 10 3 a c Jll[], providing 
nearly equipotential current distribution in the a6-plane, 
at least in the central area of the sample. Thus the con- 
tribution of a a b to the anisotropic conductivity is weak 
and can be neglected, in the limit of small pad separa- 
tion. Thus a standard 4-probe method can be used for 
p c measurement instead of the complicated multitermi- 
nal Montgomery analysis. The temperature was stabi- 
lized with an accuracy ±50 mK. The resistivity p c (Bu) 
was measured in the B± interval where g at maximum 
By drops with B\\ at least to the value 0.2. 

In Fig. 1 we present p c as a function of B 2 at different 
B± for a) irradiated and b) pristine crystals. For irradi- 
ated crystals at low Bu the resistance increases quadrati- 
cally with Bn as prescribed in Eq. (8). In contrast, for the 
pristine crystal we observed that p c increases quadrati- 
cally at high Bn , but exhibits a minimum at low fields 
that is not described by Eq. (8). This low field behav- 
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ior is caused by dynamic screening of Bu that results in 
additional dissipation at low fields, due to the combined 
effects of motion of Josephson vortices induced by Bu and 
of pancake vortices. These effects are not important at 
high fields. They do not appear in the irradiated crystal 
due to pinning of the Josephson vortices. 
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FIG. 1. Dependence of resistivity p c versus B\ for per- 
pendicular components B±, increasing sequentially with the 
step 0.1 T in irradiated (a), and pristine (b) crystals. 

In Fig. 2 we show the dependence g = a c (b) / a c (0) on 
b = B\\/ (B S B±) X / 2 for irradiated and pristine samples 
at different B±. Here B s = <&q/2its 2 . For the pristine 
crystal the values a c (B±) at -By = were determined 
by extrapolation of p c (B||) from the high field quadratic 
dependence to zero as shown in Fig. lb by the dashed 
lines. Note that for the pristine crystal all three curves 
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Fig. 2, N. Morozov 

FIG. 2. Dependence g(B n ) = a c [B n /{B s B ± ) 1/2 ]/a c (0) in 
irradiated and pristine crystals. 

coincide at high fields, demonstrating scaling of the cor- 
relation length R with a. Such a scaling is not observed 
for the irradiated crystal at B± < 1.2 T, because the 
average distance between CDs, a$ = (^o/B*) 1 / 2 , gives 
another length scale in addition to a. More importantly, 
we see that g drops with B\\ much faster for the irra- 
diated sample. The correlation length R{B±) is related 



via R w $q/B\\s to the magnetic field Bu which charac- 
terizes the scale of the drop of a c with Bu. Here Bu 
is the magnetic field at which flux in the area Rs is 
sa $ . If we define the characteristic field Bu, as given 
by: g(Bu) =0.1, then from Fig.2 at Bj_ = 0.2 T, we 
estimate for the pristine crystal R/a rj 2.5, while for ir- 
radiated one R/a « 8. As one can see from Fig. 2, at 
the high fields Bu used in our measurements, the depen- 
dence <7(-B|i) is close to exponential, and we use this to 
extrapolate <?(-B||) to higher fields. Then we determine 
G(r) in the interval r < r m « 0.5 /im by use of Eq. (|^). 
Accuracy of our experimental data and of the inverse 
Fourier-Bessel transformation is not sufficient to obtain 
G(r) at r > r m . 
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FIG. 3. The functions G{r/a) (left), and G{r/R) (right) 
as extracted from g(Bu) for irradiated (solid) and for pris- 
tine (dashed) samples. Inset: p c (B c ) at Bu = from Ref. 7. 
Arrows indicate B± values for the presented G(r) curves. 

In Fig. 3a we show the functions G(r/a) and in Fig. 3b 
the same data versus x = r/R, where the scaling length 
R(B±, B<p) is defined by Eq. 0. For the pristine crystal 
at B ± = 0.1 T and 0.2 T we obtain R(B ± ) w a. The 
curve for B± = 0.2 T is plotted by dashed line and it 
practically coincides with that at B± = 0.1 T. For the 
irradiated sample we show the curves corresponding to 
regions below, above, and at the position of the dip in 
the magnetoresistance curve, p c (B±,0) from Ref. 0] as 
shown in the inset. It is evident in Fig. 3a that the rate of 
decay of correlations with r is not a monotonic function 
of B± and is a minimum at a value corresponding to the 
dip in p c (B±). Notably all these curves quite accurately 
merge in Fig 3b, providing a universal function G(x) for 
both pristine and irradiated samples with a single scaling 
length R, although the irradiated sample is characterized 
generally by two lengths a and a$. This means that 
G(x) is determined mainly by the static effect of vortices 
on Josephson coupling 0] though it is extracted from 
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the dynamic quantity a c . In comparison with pristine 
crystals CDs simply diminish the effective concentration 
of pancakes acting on Josephson coupling by the factor 
a 2 /R 2 . 

The correlation length R as a function of B± is shown 
in Fig. 4. It exhibits a distinct maximum, R/a « 4, at 
/ w 0.35, again coinciding with the position of the dip in 
Pc(B±). At / = 0.35 the ratio of R(B±) for irradiated 
and pristine crystals is about 4 times. 
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FIG. 4. Dependence of the correlation length R as a func- 
tion of filling factor, / = B±/B<j>, of columnar defects for 
irradiated and pristine samples. Dashed lines are a guide. 

The length R± obtained from the function G(r) using 
Eq. (7) is R* 2R at all values B± studied. The same length 
determined directly from the function g(B^) at small Bn 
using Eq. (||) is « 3R. Thus we estimate the accuracy of 
extracting G(r) as s» 30 %. 

In order to determine the c-axis correlations of pan- 
cakes and to explain scaling, we note that at small / 
pancakes are positioned mainly inside CDs and hence the 
drop in the phase difference correlations is caused by in- 
terruptions in pancake arrangement inside CDs. Namely, 
we suppose that the phase difference </? n! „+i(r) is induced 
when a columnar defect in the layer n is occupied by a 
pancake, but the site inside the same CD in the layer 
n + 1 is empty (and vise versa). At larger / the un- 
pinned pancakes also contribute to the decay of corre- 
lations. Scaling for both pristine and irradiated crystals 
means that ends of vortex segments inside CDs act on the 
phase difference correlations as unpinned pancakes. The 
net concentration of interruptions and unpinned pan- 
cakes is 1/R 2 , i.e. the average length of vortex segments 
is L/s R* R 2 /a 2 . In the most aligned vortex liquid, at 
/ ~ 0.35, we estimate L/s rj 15. This is much larger 
than Lq/s rj 1/(1 — /) » 1.5 in the model of noninter- 
acting pancakes positioned randomly on CDs. Thus we 
conclude that interaction of pancakes is important for 
enhanced alignment inside CDs. 

In the hierarchy of interactions in the presence of CDs, 
both the pinning energy per pancake and the intralayer 
repulsion energy of pancakes are characterized by the 
same energy scale Eq. The scale of magnetic pair at- 



traction of pancakes in different layers is smaller by the 
factor s/A a b. The random distribution of CDs is impor- 
tant. Repulsion between vortices in the same layer leads 
to a significant increase of pancake energy inside CDs sit- 
uated near those already occupied by pancakes ]l2| |. As 
B± increases, some CDs become more favorable for filling 
by pancakes, while others remain unoccupied. Another 
important point is that favorable configurations are sim- 
ilar in all layers due to the geometry of CDs. Thus repul- 
sion of pancakes inside randomly positioned CDs in the 
same layer leads to enhancement of c-axis correlations. 
Another mechanism for pancake alignment inside CDs is 
magnetic attraction of pancakes in adjacent layers. 

In conclusion, we extract the universal phase differ- 
ence correlation function using measurements of c-axis 
resistivity as a function of the parallel component of the 
magnetic field at fixed perpendicular component. We 
estimated the correlation length of the pancake density 
correlation function along the c axis as a function of the 
filling factor / of columnar defects. It first increases with 
/, reaches a maximum at / rj 0.35, and then drops as 
pancakes start to occupy positions outside of CDs. We 
argue that enhancement of the alignment of pancakes in- 
side columnar defects is caused by interaction of pancakes 
confined inside columnar defects. 
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